Quantum state-resolved probing of strong-field-ionized xenon atoms using 
femtosecond high-order harmonic transient absorption spectroscopy 
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Femtosecond high-order harmonic transient absorption spectroscopy is used to resolve the com- 
plete \j, m) quantum state distribution of Xe + produced by optical strong-field ionization of 
Xe atoms at 800 nm. Probing at the Xe AT 4 / 5 edge yields a population distribution Pj\ m \ of 
P3/2.1/2 : Pi/2,1/2 : £3/2,3/2 = 75 ± 6 : 12 ± 3 : 13 ± 6%. The result is compared to a tunnel 
ionization calculation with the inclusion of spin-orbit coupling, revealing nonadiabatic ionization 
behavior. The sub-50-fs time resolution paves the way for table-top extreme ultraviolet absorption 
probing of ultrafast dynamics. 

PACS numbers: 42.50.Hz, 32.80.Rm, 42.65.Ky, 82.53.-k 



Studies of laser-atom interactions in the nonperturba- 
tive, strong-field regime elucidate novel phenomena such 
as above-threshold ionization [l], Q , nonsequential dou- 
ble ionization [3[, and high-order harmonic generation 
(3) B Si- While these processes are extensively stud- 
ied both experimentally and theoretically, details remain 
unknown about the \j, m) state distribution of the pho- 
toion produced by the initial photoionization step (to is 
the projection quantum number associated with the total 
angular momentum j of the hole). Moreover, experimen- 
tal tests of theoretical models for strong-field ionization 
mostly rely on measuring the ion yield as a function of 
laser peak intensity [Tj • The strong dependence of ioniza- 
tion yields on the orbital angular momentum and its di- 
rection relative to the laser polarization axis, as predicted 
by most theoretical models (e.g., the Ammosov-Delone 
Krainov (ADK) rates Q), suggests that knowledge of 
the complete |j, to) state distribution can be used as an 
additional benchmark for theory. Young et al. recently 
reported the use of synchrotron x-ray pulses to probe 
the hole-orbital alignment of Kr + photoions generated in 
the strong-field ionization of Kr [9( . The unresolved fine- 
structure transitions prevented retrieval of the complete 
\j, to) state distribution. However, the observed degree 
of alignment is reproduced by the \j, to) state distribu- 
tion obtained by tunnel ionization calculations with the 
inclusion of spin-orbit coupling 

Here we investigate the experimental and theoretical 
strong-field ionization of xenon to extract the complete 
\j, to) quantum state distribution. Femtosecond extreme 
ultraviolet (EUV) transient absorption spectroscopy is 
demonstrated with a laser-based, high-order harmonic 
probe source for the experiments; results are compared 
to tunnel ionization calculations that incorporate spin- 
orbit coupling. The resultant angular momentum distri- 
bution and hole-orbital alignment of the Xe + photoions 
are measured by probing the transition from the 4c? core 



level to the 5p valence shell. These measurements allow 
the determination of the complete \j, to) quantum state 
distribution, which is compared to theory. 

The schematic of the experimental setup is illustrated 
in Fig.[T] Briefly, the amplified output from a commercial 
Tksapphire laser system (2.4 W, 800 nm, 45 fs, 1 kHz) 
is sent to a 20 : 80 beamsplitter to produce the opti- 
cal pump and high-order harmonic generation beam, re- 
spectively. High-order harmonics in the EUV region are 
generated by focusing the laser light into a 7-cm-long, 
150-fim internal diameter capillary filled with 6.0 x 10 3 
Pa of neon [ll| . The estimated photon flux at the source 
is 10 photons per pulse for the high-order harmonic cen- 
tered at 55.4 eV. A pair of 0.2-/^m-thick Al foils is used to 
reject the residual 800-nm light and transmit the high- 
order harmonics. After reflection by a toroidal mirror, 
the high-order harmonics are refocused into a 2-mm-long 
gas cell filled with 2.7 x 10 3 Pa of Xe [13] . Scanning knife- 
edge measurements give a beam waist of 21 /jm for the 
high-order harmonics. The transmitted EUV radiation 
is spectrally dispersed in a home-built spectrometer and 
detected with a thermoelectrically-cooled CCD camera. 
A dielectric-coated mirror with a 1-mm-diamctcr intcr- 



..bearn splitter 




FIG. 1: Schematic illustration of the experimental setup. The 
light and dark lines correspond to the 800-nm pump and EUV 
probe beam paths, respectively. 
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FIG. 2: Transient absorption spectrum of Xe + acquired at a 
time delay of +500 fs, obtained from the average of 3 data sets. 
AOD denotes the change in optical density (absorbance). 



The 2 P- 
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D 5 / 2 and Pi 



/2 



-D3/2 fine-structure transi- 
tions are located at 55.4 eV and 56.1 eV, respectively. The fit 
to a sum of two Gaussian curves is depicted by the solid line. 
The inset shows the harmonic spectra obtained at —500 fs 
and +500 fs pump-probe time delay. 



FIG. 3: Time evolution of the 2 P 



3/2 



D5/2 transition for 



parallel and perpendicular relative polarizations between the 
optical pump and EUV probe pulses, obtained from the av- 
erage of 8 time scans for each relative polarization. The po- 
larization anisotropy observed at positive time delays implies 
that the 2 P%/2 state of Xe + produced by strong-field ioniza- 
tion is aligned. 



nal bore hole allows the optical pump beam to overlap 
with the EUV probe beam in a collinear geometry. A 
half-waveplate inserted into the path of the pump beam 
allows the relative polarization of the pump and probe 
beams to be varied. The polarization purity of the pump 
beam is characterized by an extinction ratio of > 200 : 1. 
The 800 nm pump pulse energy incident on the gas cell is 
0.13 mJ. Further measurements confirm that the pump 
pulse energy after the gas cell remains unchanged (by 
< 3%) upon rotation of the half-waveplate. In the pres- 
ence of the Xe gas, the pump pulse duration increases 
slightly to 49 fs and the elliptical focal spot of the pump 
beam is characterized by horizontal and vertical beam 
waists of 61 and 34 /im, respectively; these parameters 
yield a peak intensity of 8 x 10 13 W/cm 2 for the pump 
pulse. The pump-probe time delay is varied by means 
of an optical delay line in the path of the pump beam. 
Transient absorption spectra are obtained by using spec- 
tra collected at —500 fs time delay as the reference [13j; a 
negative time delay implies that the probe pulse arrives 
at the sample before the pump pulse. Error bars reported 
below correspond to 95% confidence interval limits. 

The transient absorption spectrum of Xe + acquired 
at a pump-probe time delay of +500 fs and a paral- 
lel relative polarization between pump and probe beams 
is shown in Fig. [2] [3]. The peaks at 55.4 eV and 
56.1 eV correspond to the 5p~} 2 ( 2 P 3 / 2 ) — > 4:d~} 2 ( 2 D 5 / 2 ) 



and 5pjy 2 ( P1/2) — > 4g? 3 ^ 2 ( 2 D 3 / 2 ) transitions, respec- 
tively, in agreement with literature values [HI (nlj 1 sym- 
bolizes a hole in the nl orbital with angular momentum j; 
the corresponding term symbol is given in parentheses). 
The ratio of the areas of the two fine-structure absorp- 
tion lines, defined as R = I. 
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state is followed by varying the pump-probe time delay 
while monitoring the transient absorption signal at 55.4 
eV. The resultant time traces obtained for parallel and 
perpendicular relative polarizations between pump and 
probe beams are shown in Fig. [3l Fitting the time traces 
to a convolution of a step function with a Gaussian yields 
fwhm values of 37 ± 1 fs and 39 ± 2 fs for parallel and per- 
pendicular relative polarizations, respectively. Note that 
the temporal signal corresponds to a cross-correlation of 
the Xe + population growth with the EUV pulse [l^ . 
From the polarization-dependent absorption at positive 
time delays > 50 fs, the polarization anisotropy, defined 

as r = (^3/2^5/2 _ ^3/2->5/2)/(^3/2->5/2 + 2 ^3/2— 5/2)' is 

found to be r = 0.07±0.01. The observed anisotropy im- 
plies the existence of hole-orbital alignment in the Xe + 
2 P 3 / 2 state produced by strong-field ionization. Since 
the hole orbital is directed along the polarization axis 
of the pump beam (which also defines the quantization 
axis), the hole population of the m — ±1/2 sublevels is 
expected to be greater than that of the m = ±3/2 sub- 
levels, as observed experimentally. To verify the reliabil- 
ity of the measured polarization anisotropy, a separate 
set of measurements is performed for the 2 Pi/2 — > 2 D 3 / 2 
transition, for which no anisotropy is observed. This re- 
sult agrees with the fact that alignment cannot exist in 
a 2 Pi/ 2 state (Tt | . 



The experimental results for the ratio of the fine- 
structure absorption R and the polarization anisotropy r 
can be used to extract the complete quantum state distri- 
bution. The description of the EUV-probe step follows 
Ref. Ell. Let pjA m \ denote the probability of finding 



is found to Xe + with a hole in either the 5»j m or the bpj 



be R = 6.5 ± 1.1. The time-evolution of the Xe+ 2 P 



3/2 



bital. Making the dipole approximation and employing 
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standard angular momentum algebra [17} , we obtain 

_ 1 P3/2A/2 — P3/2,3/2 



and 



R 



10 P3/2.1/2 + P3/2.3/2 

_ 3p 3 / 2 ,l/2 + 2p 3 / 2 ,3/2 ^ 
5^1/2,1/2 



(1) 



TABLE I: Comparison of the complete \j, m) quantum state 
distribution obtained from experiment and theory for Xe + 



(2) 



where £ = \( 2 D 5/2 \\ d \\ 2 P 3/2 )\ /\( 2 D 3/2 \\ d \\ *P 1/2 )\ 2 
and d is the electric dipole operator. A multiconfigura- 
tion Dirac-Fock calculation performed with the program 
package GRASP2 [l| gives £ = 1.6. From Eqs. O 
and ([2]), and the experimental values for R and r, the 
complete quantum state distribution of Xe + generated 
by strong-field ionization can be extracted. The results 
are summarized in Table HI We note that spectroscopic 
probing of the photoion by transient absorption allows 
direct retrieval of its quantum state distribution, which 
cannot be obtained from other experimental approaches 
to the study of strong-field ionization, such as energy- 
and angular-resolved photoelectron spectroscopy [3] . 

Tunnel ionization calculations with the inclusion of 
spin-orbit coupling are performed to model the exper- 
imental results described above. To calculate the pro- 
duction rates for the various ionization chann els, the ef- 
fective one-electron model described in Ref. [lfj is em- 
ployed. The model treats strong-field ionization within 
the tunneling picture and includes the effect of spin- 
orbit interaction. The approach makes use of a flexible 
finite-element basis set and determines ionization rates 
in this square-integrable basis using a complex absorb- 
ing potential [23]. Parameters that are used in the cal- 
culations are 0.13 mJ for the pulse energy, 49 fs fwhm 
for the pulse duration, and 45 /im for the cylindrically- 
symmetric beam waist. The \j, m) populations are cal- 
culated after the laser pulse by numerical integration 
of rate equations [Io| and are normalized such that 
P3/2.1/2 + Pi/2,1/2 + P3/2.3/2 = 1- The calculated frac- 
tional populations on the laser axis are shown in Table [J 
These calculations also allow us to verify that effects due 
to spatial averaging by the probe beam are smaller than 
the experimental errors for the Pj\ m \ distribution. 

Within the given uncertainties, good agreement for the 
P3/2.1/2 : Pi/2,1/2 ratio is obtained between experiment 
and theory. This suggests that spin-orbit interaction is 
adequately treated in the tunnel ionization model em- 
ployed in the calculations. Moreover, the calculations re- 
veal that there is very little mixing between the 5^3/2 and 
5pi/ 2 valence orbitals even at the saturation intensity for 
Xe + production. This implies that the laser field is not 
sufficiently strong to uncouple the spin-orbit interaction. 

Furthermore, the tunneling calculation predicts that 
P3/2,3/2 *C P3/2.1/2) which supports the experimental ob- 
servation that the Xe + 2 P 3 / 2 state produced by strong- 
field ionization is aligned. However, it is notable that 
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the measured ^3/2,3/2 : P3/2,i/2 ratio of 0.17 ±0.09 is sig- 
nificantly larger than that predicted by the calculation, 
which gives a value of 0.04 for this ratio. The discrep- 
ancy between experiment and theory can be attributed 
to the partial breakdown of the adiabatic (quasi-static) 
approximation employed in the tunnel ionization model. 

The adiabatic approximation requires that the elec- 
trons within the atomic potential and those undergoing 
tunnel ionization respond instantaneously to the laser 
field. This approximation is analogous to the Born- 
Oppenheimer approximation frequently invoked in the 
study of molecular dynamics, whereby electrons are as- 
sumed to respond instantaneously to the electric field 
exerted by the nuclei. The ratio of the tunneling time to 
the period of the laser field is given by the Keldysh adi- 
abaticity parameter 7 plj ], implying that the adiabatic 
approximation is valid only when 7 <C 1. Given that the 
experimental conditions yield 7 = y/ I P /2U P ~ 1.1 (I p is 
the atomic ionization potential and U p is the pondero- 
motive potential), such a quasi-static approximation is 
no longer wholly valid, i.e., ionization becomes nonadia- 
batic with respect to the laser field. Under such circum- 
stances, a nonperturbative multiphoton Floquet treat- 
ment of strong-field ionization becomes necessary [22j ]. 
Indeed, previous investigations of Xe strong-field ioniza- 
tion revealed discrete peaks atop above-threshold ion- 
ization features in the photoelectron spectra, suggesting 
that ionization occurs mainly via a multiphoton path- 
way [23l . I24] ] ; the discrete peaks originate from resonance- 
enhanced multiphoton ionization that occurs when the 
field-dressed ground state crosses the ac Stark-shifted 
Rydberg states [25]. On the sub-cycle timescale, nona- 
diabatic transitions can occur between the field-dressed 
ground state and the Rydberg levels while varying the 
instantaneous phase of the laser field, akin to transi- 
tions between adiabatic potential energy surfaces in the 
non-Born-Oppenheimer regime of molecular dynamics. 
Note that nonadiabatic electron dynamics have previ- 
ously been observed in molecular strong-field ionization, 
in which the spatial derealization of electrons leads to 
the breakdown of the adiabatic approximation [26j . 

In the present work, the importance of nonadiabatic ef- 
fects in atomic strong-field ionization is substantiated by 
the results of two spin-orbit-free calculations — a tunnel- 
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ing calculation and a multiphoton Floquet-type calcula- 
tion 27J — performed at an intensity of 8 x 10 13 W/cm 2 . 
The Floquet calculation predicts that the ratio of the 
rrii = il : m i = ionization rates is ~ 2x larger than 
that predicted by the tunnel ionization model. There- 
fore it is reasonable that the experimentally-measured 
/ 9 3/2,3/2 : Pz/2,1/2 ratio is larger than that predicted by 
the tunneling calculation with spin-orbit coupling. A 
valid theoretical comparison with the experimental re- 
sult would require a Floquet calculation that incorporates 
spin-orbit coupling. 

Finally, we note that while recent advances in high- 
order harmonic generation have already resulted in the 
extension of ultrafast spectroscopy into the EUV do- 
main 111 3, 3(3, 1H, most of the experiments reported 
to date have focused on photoelectron spectroscopy due 
to the set of discrete harmonics that presents itself as 
an attractive photoionization source. By probing transi- 
tions from the core level to unoccupied valence levels, pi- 
cosecond time-resolved x-ray absorption near-edge spec- 
troscopy based on laser-produced plasma [13] and syn- 
chrotron 33j, l34j sources have been shown to be highly 
sensitive to element-specific electronic structure changes 
accompanying photophysical and photochemical trans- 
formation. The work here demonstrates the feasibility of 
performing transient absorption spectroscopy with sub- 
50-fs time resolution using a laser-based, high-order har- 
monic source for core-level probing. Work on extending 
this technique to the study of ultrafast molecular dynam- 
ics is currently in progress. 
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